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A consistent BPS formalism to study the existence of topological axially symmetric vortices in 
generalized versions of the Born-Infeld-Higgs electrodynamics is implemented. Such a generalization 
modifies the field dynamics via introduction of three non-negative functions depending only in 
the Higgs field, namely, G(|())|), ic(|0|) and H(|<(>|). A set of first-order differential equations is 
attained when these functions satisfy a constraint related to the Ampere law. Such a constraint 
allows to minimize the system energy in such way that it becomes proportional to the magnetic 
flux. Our results provides an enhancement of topological vortex solutions in Born-Infeld-Higgs 
electrodynamics. Finally, we analyze a set of models such that a generalized version of Maxwell- 
Higgs electrodynamics is recovered in a certain limit of the theory. 

PACS numbers: ll.lO.Kk, ll.lO.Lm, 11.27.-|-d 


I. INTRODUCTION 

The well-known Born-Infeld electrodynamics was orig¬ 
inally introduced to remove the divergence of electron’s 
self-energy in classical electrodynamics by introducing 
a square-root form of the Lagrangian density replacing 
the standard Maxwell Lagrangian [I|. In this way the 
field strength tensor remains bounded everywhere and 
the energy associated to a point-like charge becomes fi¬ 
nite. This theory is a distinguished member of the family 
of nonlinear electrodynamics since it enjoys three prop¬ 
erties: i) Maxwell weak-field limit, ii) electric-magnetic 
duality (2|, and iii) absence of shock waves and birefrin¬ 
gence phenomena concerning propagation of waves, be¬ 
longing to the class of theories called “completely excep¬ 
tional” 0 . Applications of Born-Infeld electrodynamics 
within gravitation and cosmology have been considered 
for many years Q. This model is moreover deemed of a 
special attention since it appears in the low-energy limit 
of string/D-Branes physics 

On the other hand, the study of magnetic vortices 
gained great interest since Abrikosov’s description for 
Type-II superconductors which arise naturally from 
the non-relativistic limit of Ginzburg-Landau (GL) the¬ 
ory 0 - In field theory, stable vortex configurations came 
up with the seminal work by Nielsen and Olesen M 
whose study of the Maxwell-Higgs (MH) model shows 
that electrically neutral vortex solutions correspond to 
the ones obtained by Abrikosov. Lately it was verified 
the existence of electrically charged vortex solutions in 
Chern-Simons-Higgs (GSH) [ij, and Maxwell-Ghern- 
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Simons-Higgs (MCSH) models. In all these cases, 

the presence of the Higgs fields is essential for the exis¬ 
tence of vortex-type solutions. 

Recently, it has been intensively studied the existence 
of topological defects in generalized or new effective field 
theories. For example one can introduce noncanonical 
kinetic terms [l^ . Ilq . in order to circumvent the con¬ 
straints of Derrick’s theore m 11^ and obtain topological 
defect solutions (see e.g. [l^ for a more detailed ac¬ 
count on soliton-like solutions in field theory). Other 
models are defined by introducing generalizing functions 
on standard field models M- In some cases these gener¬ 
alized models provide self-dual analytical solutions which 
certainly enriches our understanding of the field [13, [HI . 
Moreover, this procedure allows to control properties of 
the topological defect, such as its width or energy den¬ 
sity, providing valuable models for the analysis of several 
physical problems. In the literature there are many in¬ 
teresting applications of these new solutions within sev¬ 
eral different scenarios, specially involving the acceler¬ 
ated inflationary phase of the universe [22| via the so- 
called k-essence models 23|, strong gravitational waves 
[23| . t achy on matter [25 1, dark matter [l^l, and other 
topics [27f . 

Among generalized models the simplest ones are 
those generalizing the Maxwell-Higgs model [l^l , Chern- 
Simons-Higgs model [131 and Maxwell-Chern-Simons- 
Higgs model [131. Based on earlier work on vortices in 
Born-Infeld-Higgs models [HI, in Ref. [l3| a general¬ 
ization of Born-Infeld-Maxwell-Higgs (BIMH) model was 
constructed within the context of generalized dynamics, 
but self-dual or BPS vortices were not found. The main 
aim of the present manuscript is to show the existence of 
self-dual topological BPS vortices in a generalized BIMH 
electrodynamics, and study their properties. 
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II. GENERALIZED BORN-INFELD VORTICES 

The Lagrangian density of our (2 + l)-dimensional the¬ 
ory is written as 

c = p^{i-n) + w (|,^|) - w (|<^|), ( 1 ) 

with the definitions 

n = + ( 2 ) 

wm = ^^[i-vm], ( 3 ) 

where = d^A^, — d^A^j, is the field strength tensor 
of the vector potential A^, while the covariant deriva¬ 
tive realizing the coupling between the gauge and Higgs 
fields is given by = 9^(/) — leA/^cj). The positive 
functions Cdf/il) and w {\(j)\) are the generalizing func¬ 
tions in the kinetic sector. The generalized potential 
IT (|(/)|), a nonnegative function, inherits its structure 
from the function V(|(/)|), which is restricted by the con¬ 
dition 0 < T(|(/)|) < 1, so W[(j)) > 0. The Born-Infeld 
parameter, j3, provides modified dynamics for both scalar 
and gauge fields further enriching the family of possible 
models. 

From the action 0 the gauge field equations of motion 
read 

d. ( 4 ) 

where — 2eA>^ plays the role 

of a current. 

At static regime, Eq. dU provides the Gauss law, 

dj = ‘2‘We^Ao , ( 5 ) 

which is saturated by temporal gauge: Aq = 0. In this 
way we see that the model, at static regime and in tem¬ 
poral gauge, describes electrically neutral magnetic con¬ 
figurations. Under these conditions, from Eq. @ the 
Ampere law reads 


ekj dj 



ewJk = 0 , 


and Higgs’s field equation becomes 


( 6 ) 


The energy-momentum tensor of the model is given by 
= -^G M) + w (|<^|) {D^4>r D,cj3 

+w{\cl)\){D^(j))* D^cj)- (9) 

In this work we are interested in searching for electri¬ 
cally neutral magnetic vortices and, more specifically, we 
will study such solutions at static regime and in tempo¬ 
ral gauge. As it is largely known in literature, the axially 
symmetric vortex ansatz works fine to find such solutions, 
namely, 

4> = vg{r)F-\ = ( 10 ) 

er 

where n is an integer number and a(r) and g{r) are reg¬ 
ular functions that satisfy the following boundary condi¬ 
tions 


5 ( 0 ) = 0, 5 ( 00 ) = I (11) 

a ( 0 ) = n , a ( 00 ) = 0 ( 12 ) 

Using this ansatz the magnetic field is written as 

B = --. (13) 

er 

with the short-hand notation a' = da/dr. 

For the ansatz (nni) the Ampere law m is expressed as 

=-2ev^w^, (14) 








^ ,, g' a? 1 dw 

+ 7 - 723 +^ 


2w dg 


w 


I dG 1 dW 

4wv'^ TZ dg 2wv‘^ dg ’ 

where TZ is given by Eq.®. 



(15) 


A. The BPS formalism 

The energy of the vortex is given by the integration of 
the Too component of ® which, in static regime and in 
the gauge Ag = 0 , is given by 


0 = w {DjDjip) + (djw) Djcf) 


dw 2 dG dW 


where in the last two equations TZ reads 




1/2 


C^) Too = {TZ-V) + wv'^ {g'f+ wv'^ , (16) 

and will be nonnegative whenever the condition TZ > V 
is satisfied. The total energy reads 

E = 2n 

( 8 ) 


Jdrr /3^ {TZ — V) + wv^ {g')^ + wv^ (~) 


TZ = 


( 17 ) 
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where the fields were expressed in terms of the ansatz 
(uni). We now use the Bogomol’nyi trick [32j to rewrite 
it as 



where we have introduced the function F, which is, in 
principle, arbitrary but nonnegative, to be determined 
later in order to obtain solutions with well defined energy. 
Using the definition ([8]) in the third row, we can rewrite 
(fT51) as 

E = 2Trjdrrl±B^y/WG±2wv^'^g' (19) 


where 

with the function F given by Eq. (Uni- We can now 
further constrain the set of functions G, w and V in order 
to attain a true lower bound for the energy by selecting 
functions satisfying 

n (0) = 1, U (oo) = finite (25) 

Then, by considering the boundary conditions given by 
Eqs. (EH) and (HU, the energy reads 

E = 2-kv'^ \n\ (26) 



This clearly shows that the energy possess a lower bound 



We observe that, by imposing the expression in the third 
row to be null, this allows to determine the function F 
in terms of V and R, namely 


F = ^{n-vf + ^{i-v^), ( 20 ) 


which shows that E is a nonnegative function because 
0 < U < 1. Let us point out that the function F is 
defined without considering the self-dual equations or the 
BPS limit. 

Now by considering condition (1201) and the expression 
m for the magnetic field, the energy (ITOl) reads 



E > 2t:v^ \n \, (27) 

whenever the functions G, w and V chosen provide a 
function il satisfying Eq. (P5|) . Such a lower bound is 
saturated when the fields satisfy the BPS or self-dual 
equations 



(28) 



1-U2 

G 


(29) 


This is a set of first-order equations that satisfy automat¬ 
ically the second-order equations da and m, as can be 
immediately seen by derivation of the former. This is so 
because the Euler-Lagrangian equations only imply that 
a static BPS solution will be a stationary point of the 
energy. In Ea. l^H)) we have used (1^ to compute F in 
the BPS limit, which provides 


1 - 


21/2 ’ 


(30) 


a nonnegative function due to 0 < U < 1. Similarly, the 
nonnegative function (r) is given by 


We can now transform the first term in a total deriva¬ 
tive by setting 

=-2wgg', (22) 

This way, the energy dUD is written as 



flbp. = /3VG(1-U2). (31) 

By using the BPS equations (l28l) and (l29l) . Ampere’s 
law da can be written as 

-^\//32G (1 - U2) = -2ev^wg. (32) 

dg 

This relation allows to determine one of the generalizing 
functions when the other two are given, for example, we 
can compute w if we provide the functions G and V. 
Here it is worthwhile to notice that Ea. (|5^ is exactly 
the condition (1^^ in the BPS limit. 
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To conclude this section, the BPS energy density of 
the model, which appears in 

^bps — ‘Ztt J^dvr £f)ps^ 

is given by 

Sbps = . (34) 

and will be positive definite whenever the functions 0 < 
( 5 ) < 1 and w [g) > 0. 

III. A FAMILY OF MODELS 


will consider, in each case, the asymptotic behavior of 
the functions g{r) and a(r) compatible with the bound¬ 
ary conditions that make the energy finite and positive, 
and next solve the BPS equations. On the other hand we 
note that f3 is not a constant characterizing the solutions 
but rather a parameter determining a particular model 
within the family defined by the corresponding term in 
the action (HD. In some of the following numerical cases 
we shall treat nevertheless /3 as a free parameter for the 
computations, which means that in those cases we will be 
comparing the behavior of the solutions corresponding to 
different models of the family of generalized Born-Infeld 
Lagrangians. In order to perform the numerical analysis, 
without loss of generality, we set e = 1 = v. 


In this section we shall focus on the special case 


B. Asymptotic behavior for models 


VM) 


'1- 


2 um 


(35) 


The \(j)\ -models are described by the function U [g) 
given by 


since this choice, in the limit /3 —>■ 00 , allows to obtain the 
generalized Maxwell-Higgs model from the Lagrangian 
density ([T]): 

c = + wm) - U (|</.|). (36) 


U{g) = ^{l-g^)\ (41) 

and the function G {g) whose behavior when r —>■ 0 is 

G {g) = ao + a 2 g^ + ■■■, (42) 


With this choice the BPS equations read 



The condition (l32l) reads 

-^V2t/G = -2ev^wg, (39) 

dg 

and the BPS energy density is 

ebps=2U +2^^"w(y) • (40) 

Therefore, the generalized models can be defined by 
choosing G {g) and U (g) functions which, via the con¬ 
straint dMD, allow to find the remaining function w (g). 
These three functions must be nonnegative for positive 
definiteness of the energy density. In the next sections 
we shall choose some models satisfying the constraint 
(1^51) and, therefore, their BPS solutions will saturate the 
bound 071) . 


A. Some choices for the potential 


and when r —> 00 reads 

G (g) = (1 - 5) + ^2°°^ (1 - gf + (43) 


where ao, a 2 ,..and aQ°°\ ... are some constants. 

By introducing the above information into the BPS 
equations (1^ and (IMl) , we can compute the behavior of 
field profiles when r —)> 0: 


g (r) = Cnr" 


a (r) = n — 



1/2 G 


V/3"-V 4(ao)'^^ 

yG ^2 

/3^-iy 2{a,f^ 

(2q!o -b 02 ) /?^ - 02 G^r^"+^ 
4(n-bl) 


(44) 

(45) 

-b ... 


where G„ is a set of constants. 

Similarly, we calculate the behavior of the profiles at 
infinity: 


1 - 5 (»■) 
a{r) 



(46) 

(47) 


Next we shall consider two classes of models charac¬ 
terized by the form of the “potential” U{g). First we 


These expansions are fully consistent with the assumed 
boundary conditions (ED and (|12l) for the BPS solutions. 
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C. Asymptotic behavior for \(j)\ models 


In this case, the function U (g) is given by 


(48) 


We consider the behavior of function G, which at origin 
takes the form 


G {g) = ^ + 70 + 725^ + 
9 


(49) 


and at infinity is 

G {g) = 7o°°^ + 7 !°°^ (1 - 5 ) + 72°°^ (1 - gf + - (50) 

The behavior of the profiles at r —0 is 

(j3j,3n+2 


9ir)r^G„r^- 
a(r) ^ n — 


4 ( 7 - 2 )^^^ (n+ 1)^ 

(j2^2n+2 


2 ( 7 - 2 )'/" (n + 1) 

, (27-2 +7o)/ 3^-7-2 ^4„4n+2 , 

I Q /o ^ n' 

4/32 ( 7 - 2 ) / ( 2 n+ 1 ) 

while the asymptotic behavior for r —>■ 00 is 

, -1/4' 


(51) 

(52) 


1 — g (r) ^ r 0 xp 
a (r) ^ gxp 


-r\f2 ( 70 “^) 


(53) 

(54) 


and, again, these expansions are consistent with the prob¬ 
lem under consideration. 



Figure 1. Higgs field g(r), for n — 1 and dashed-dotted lines 
for /3 = 1.5. 



Figure 2. Gauge field profile a(r), for n = 1 and dashed- 
dotted lines for P = 1.5. 


D. Discussion of results 

Once the boundary conditions are fixed, we have per¬ 
formed numerical solutions of the BPS equations dSZD 
and ((551) by using routines of Maple 16.2. The first nu¬ 
merical results are obtained by considering fixed values 
of /3(= 1.05), and comparing the standard MH, CSH and 
BIMH models with our |())|"'-BIMH and |(/)|®-BIMH mod¬ 
els. These results are shown in Figs. (U [5115] and S) The 
second numerical analysis was performed by fixing n = 1 
and varying the values of /3(= 1.05,1.25, 2.00, 00 ), with 
the resulting profiles depicted in Figs. [51111 [3and|5]for the 
|(/)|"'-BIMH and |0|®-BIMH models studied in this work. 
In both scenarios, we have depicted the field profiles g{r), 
a(r), the magnetic field B{r) and the BPS energy den¬ 
sity ebps{r) corresponding to the different models under 
comparison. 

To further clarify the plots, we note that the first 
model is defined by 

(55) 


and represents the standard Born-Infeld-Maxwell-Higgs 
model (dashed-dotted red lines in Figs. HHH). 

The second one (dashed-dotted green lines in Figs. [B- 
Sl) is given by the following functions 

G (g) = exp {2g^) , w {g) = g'^ exp {2g‘^) (56) 

The (/)®-model (dashed-dotted orange lines in Figs. (BSl) 
is defined by the functions 

^ ( 5 ) = q^2 ^ ( 5 ) = ^ (1 + ff") ■ (57) 

For completeness, we also depict the profiles of the stan¬ 
dard MH (solid black line) and CSH models (solid blue 
line). 

In general we see that the introduction of a finite value 
for P has a non-trivial impact on the profiles of a(r) and 
g{r). This follows from the comparison between the stan¬ 
dard (j)'^ BIMH model in Eq. ([55]) (red dashed curve, 
corresponding to /3 = 1.5) and the standard MH system 
(solid black curve), with the former vortex being ticker 


G = 1, li) = 1, 
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Figure 3. Magnetic field profiles B{r), for n = 1 and dashed 
lines for /3 = 1.5. 



r 


Figure 4. BPS energy density £^p^{r), for n = 1 and dashed 
lines for /3 = 1.5. 


than the latter. We also see that the impact of changing 
the G and w functions through the new |(/)|'^ and 
BIHM models introduced in this work is to made the 
vortex even more thicker (green and orange curves, cor¬ 
responding to models (l56l) and (l57l) . respectively). This 
is also reflected in the physical magnitudes characteriz¬ 
ing the vortex, as both the magnetic field and the energy 
density profiles (see Figs. |3]and|H respectively) undergo 
large modifications as compared to their standard coun¬ 
terparts. In general this means that, at fixed /3, one can 
control thickness and physical magnitudes of the vortex 
by introduction of suitable G and w functions. 

Hereafter, we depict the profiles for the second and 
third models by fixing n = 1 and some values of (3. From 
Figs. [S]and|n]we see that for |(/)|"'-BIMH model the thick¬ 
ness of the vortex increases as (3 decreases, i.e., when 
the nonlinear effects of the Born-Infeld contribution grow 
stronger, while for the |0|®-BIMH model the new effects 
play a very little role, leaving almost unmodified the vor¬ 
tex profile. For the |(()|^-BIMH model this implies large 
modifications on the magnetic field and energy density 




r 


Figure 5. Higgs field profile g{r), for n = 1. Upper figure 
represents the |(/)|^-model defined by Eq. (I56|l and the bottom 
figure represents the |(?!>|®-model defined in Eq. 




0123456789 

r 


Figure 6. Gauge field profile a{r), for n = 1. Left figure 
represents the |</)|'*-model defined by Eq. (I56|l and right-figure 
represents the |(;fi|®-model defined in Eq. (1571) . 


profiles, since their maximum at r = 0 grows quickly with 
1//3. On the other hand, as one could have expected, the 
tin^ modifications on the vortex profile with /3 in the 
\(f>\ -BIMH model also have little effect on the magnetic 
field and energy density profiles. For this model these 
profiles have a different behavior as in the |^|"'’-BIMH 
model, since their maxima are not attained at r = 0, but 
rather at a finite distance, a feature that holds for any 
value of f3. 

This analysis shows that the modified-BIHM models 
through (3 corrections do not change the qualitative fea¬ 
tures of the physical magnitudes characterizing the vor¬ 
tex, but are able to introduce quantitative modifications. 
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Figure 7. Magnetic field B{r), for n = 1. Left figure rep¬ 
resents the |</)|'*-model defined by Eq. (1561) and right figure 
represents the |0|®-model defined in Eq. (15711 . 



Figure 8. BPS energy density Ehpsir), for n = 1. Left figure 
represents the |0|'‘-model defined by Eq. (1561) and right figure 
represents the |0|®-model defined in Eq. (1571) . 


which can become large, as in the |0|^-BIMH model. 


IV. CONCLUSIONS 


vortices has a lower bound 

E > 2Trv^ \n\ , (58) 

which is saturated by the self-dual or BPS topologi¬ 
cal solutions. In the numerical analysis we have em¬ 
ployed two classes of models characterized by the po¬ 
tential term, namely, and models, and de¬ 

picted the corresponding results for the field profiles 
and the physical magnitudes characterizing the vor¬ 
tices. Such results have been compared to those of the 
standard Maxwell-Higgs, Chern-Simons-Higgs and Born- 
Infeld-Maxwell-Higgs models. 

As observed in other cases of Born-Infeld-type modifi¬ 
cations in the literature, the introduction of finite values 
for the Born-Infeld (3 has a non-trivial impact on the 
field profiles of the vortices, with the result that the cor¬ 
responding physical properties can be controlled by ade¬ 
quate combination of Born-Infeld modification and uj{g) 
and G{g) functions. When we vary /3, however, the size 
of the variation of the vortex properties largely depends 
on the model chosen, with the |(()|^ one showing impor¬ 
tant variation, while the |(^|®-one is almost insensitive to 
changes in /3. Since topological defects find applications 
to many context of modern physics as useful tools for the 
modelling of different kinds of systems, to be able to mod¬ 
ify the physical properties of vortex solutions is a strong 
motivation in favour of consideration of this kind of mod¬ 
els. Finally let us mention that the parameter P can not 
be made arbitrarily small. Our numerical analysis shows 
that for all |(()|^-models the solutions are obtained when 
P > 1. In the case of the |(/)|®-models it was observed 
that when /3 > 1 the numeric computations are always 
valid. The presence of a critical minimum value. Pc, be¬ 
low which numerical computations break down and no 
solution can be attained, seems to be a quite general phe¬ 
nomenon occurring in Born-Infeld type modifications, as 
found in other investigations in the literature [Ig, . 

In those cases, around Pc the physical magnitudes char¬ 
acterizing the topological defect change abruptly as p is 
slowly varied, as happens in our case. Though some re¬ 
search has been performed about the implications of this 
feature, this issue remains unsolved. To conclude, we 
point out that the results presented here could be gener¬ 
alized to include non-symmetric BPS fields. 


In this work we have studied a family of generalized 
Born-Infeld theories with a free parameter, /3, and three 
generalizing functions which are nonnegative. These gen¬ 
eralizing functions are constrained by the condition (1321) 
which is the Ampere law of the model. We have worked 
out the theory and obtained BPS solutions of vortex- 
type using Bogomol’nyi trick and determined the phys¬ 
ical properties of the solutions in terms of the magnetic 
flux and energy density. It was shown that whenever the 
conditions (1^ are satisfied, the energy of the topological 
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